home *** CD-ROM | disk | FTP | other *** search
/ TeX 1995 July / TeX CD-ROM July 1995 (Disc 1)(Walnut Creek)(1995).ISO / web / fweb / fweb-1.40 / demos / Newton.f < prev    next >
Text File  |  1993-10-29  |  2KB  |  89 lines
  1.  
  2. C FTANGLE v1.30, created with UNIX on "Tuesday, May 11, 1993 at 10:55." 
  3. C COMMAND LINE: "ftangle Newton.web -z"
  4. C RUN TIME: "Thursday, June 10, 1993 at 19:55."
  5. C WEB FILE:    "Newton.web"
  6. C CHANGE FILE: (none)
  7. C* 7: * 
  8. *line 162 "Newton.web"
  9.       program cosine
  10.       implicit undefined(a-z)
  11. C* 8: * 
  12. *line 179 "Newton.web"
  13.       real*8 x_0,x,delta_x
  14. C* :8 * 
  15. C* 9: * 
  16. *line 185 "Newton.web"
  17.       integer k,limit
  18. C* :9 * 
  19. C* 14: * 
  20. *line 250 "Newton.web"
  21.       real*8 delta_x_max,delta_x_prev
  22. C* :14 * 
  23. *line 165 "Newton.web"
  24.       
  25. C* 10: * 
  26. *line 200 "Newton.web"
  27.       x_0=1.2d0
  28.       delta_x=0.000001d0
  29.       limit=10
  30. C more than safe 
  31. C* :10 * 
  32. C* 15: * 
  33. *line 255 "Newton.web"
  34.       
  35.       delta_x_prev=0.0d0
  36.       delta_x_max=0.5d0
  37. C* :15 * 
  38. *line 166 "Newton.web"
  39.       
  40. C* 11: * 
  41. *line 210 "Newton.web"
  42.       k=1
  43.       do while((delta_x.GT.0.0d0).AND.(k.LE.limit))
  44. C* 12: * 
  45. *line 231 "Newton.web"
  46.       delta_x=-(cos(x_0))/(-sin(x_0))
  47. C* :12 * 
  48. *line 213 "Newton.web"
  49.       
  50. C* 16: * 
  51. *line 264 "Newton.web"
  52.       if(delta_x.GT.delta_x_max)then
  53.       delta_x=delta_x_max
  54.       else if(delta_x.LT.(-delta_x_max))then
  55.       delta_x=-delta_x_max
  56.       end if
  57.       
  58. C* :16 * 
  59. C* 17: * 
  60. *line 281 "Newton.web"
  61.       if((delta_x*delta_x_prev).LT.0.0d0)then
  62.       write(*,*)'Oscillating',k
  63.       end if
  64. C* :17 * 
  65. *line 214 "Newton.web"
  66.       
  67. C* 13: * 
  68. *line 237 "Newton.web"
  69.       x=delta_x+x_0
  70.       x_0=x
  71.       delta_x_prev=delta_x
  72. C* :13 * 
  73. *line 215 "Newton.web"
  74.       
  75.       k=k+1
  76.       end do
  77. C* :11 * 
  78. *line 167 "Newton.web"
  79.       
  80. C* 18: * 
  81. *line 290 "Newton.web"
  82.       write(*,*)'The solution to cos(x)=0 is ',x
  83. C* :18 * 
  84. *line 168 "Newton.web"
  85.       
  86.       stop
  87.       end
  88. C* :7 * 
  89.